Analysis of GLDS-46 from NASA GeneLab
This R markdown file was auto-generated by the iDEP website Using iDEP 0.91, originally by Steven Xijin.Ge@sdstate.edu
Ge SX, Son EW, Yao R: iDEP: an integrated web application for differential expression and pathway analysis of RNA-Seq data. BMC Bioinformatics 2018, 19(1):534. PMID:30567491
First we set up the working directory to where the files are saved.
setwd('~/Documents/HTML_R/GLDS46')
R packages and iDEP core Functions. Users can also download the iDEP_core_functions.R file. Many R packages needs to be installed first. This may take hours. Each of these packages took years to develop.So be a patient thief. Sometimes dependencies needs to be installed manually. If you are using an older version of R, and having trouble with package installation, try un-install the current version of R, delete all folders and files (C:/Program Files/R/R-3.4.3), and reinstall from scratch.
if(file.exists('iDEP_core_functions.R'))
source('iDEP_core_functions.R') else
source('https://raw.githubusercontent.com/iDEP-SDSU/idep/master/shinyapps/idep/iDEP_core_functions.R')
We are using the downloaded gene expression file where gene IDs has been converted to Ensembl gene IDs. This is because the ID conversion database is too large to download. You can use your original file if your file uses Ensembl ID, or you do not want to use the pathway files available in iDEP (or it is not available).
inputFile <- 'GLDS46_Expression.csv'
sampleInfoFile <- 'GLDS46_Sampleinfo.csv'
gldsMetadataFile <- 'GLDS46_Metadata.csv'
geneInfoFile <- 'Arabidopsis_thaliana__athaliana_eg_gene_GeneInfo.csv' #Gene symbols, location etc.
geneSetFile <- 'Arabidopsis_thaliana__athaliana_eg_gene.db' # pathway database in SQL; can be GMT format
STRING10_speciesFile <- 'https://raw.githubusercontent.com/iDEP-SDSU/idep/master/shinyapps/idep/STRING10_species.csv'
Parameters for reading data
input_missingValue <- 'geneMedian' #Missing values imputation method
input_dataFileFormat <- 1 #1- read counts, 2 FKPM/RPKM or DNA microarray
input_minCounts <- 0.5 #Min counts
input_NminSamples <- 1 #Minimum number of samples
input_countsLogStart <- 4 #Pseudo count for log CPM
input_CountsTransform <- 1 #Methods for data transformation of counts. 1-EdgeR's logCPM 2-VST, 3-rlog
readMetadata.out <- readMetadata(gldsMetadataFile)
library(knitr) # install if needed. for showing tables with kable
library(kableExtra)
kable( readMetadata.out ) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| ATM1_HZE_15hpostIR_5do_Rep1 | ATM1_HZE_15hpostIR_5do_Rep2 | ATM1_HZE_12hpostIR_5do_Rep1 | ATM1_HZE_12hpostIR_5do_Rep2 | ATM1_HZE_24hpostIR_6do_Rep1 | ATM1_HZE_24hpostIR_6do_Rep2 | ATM1_HZE_6hpostIR_5do_Rep1 | ATM1_HZE_6hpostIR_5do_Rep2 | ATM1_noIR_CTRL_5do_Rep1 | ATM1_noIR_CTRL_5do_Rep2 | ATM1_noIR_CTRL_6do_Rep1 | ATM1_noIR_CTRL_6do_Rep2 | ATM1_YR_15hpostIR_5do_Rep1 | ATM1_YR_15hpostIR_5do_Rep2 | ATM1_YR_12hpostIR_5do_Rep1 | ATM1_YR_12hpostIR_5do_Rep2 | ATM1_YR_24hpostIR_6do_Rep1 | ATM1_YR_24hpostIR_6do_Rep2 | ATM1_YR_6hpostIR_5do_Rep1 | ATM1_YR_6hpostIR_5do_Rep2 | WT_HZE_15hpostIR_5do_Rep1 | WT_HZE_15hpostIR_5do_Rep2 | WT_HZE_12hpostIR_5do_Rep1 | WT_HZE_12hpostIR_5do_Rep2 | WT_HZE_24hpostIR_6do_Rep1 | WT_HZE_24hpostIR_6do_Rep2 | WT_HZE_3hpostIR_5do_Rep1 | WT_HZE_3hpostIR_5do_Rep2 | WT_HZE_6hpostIR_5do_Rep1 | WT_HZE_6hpostIR_5do_Rep2 | WT_noIR_CTRL_5do_Rep1 | WT_noIR_CTRL_5do_Rep2 | WT_noIR_CTRL_6do_Rep1 | WT_noIR_CTRL_6do_Rep2 | WT_YR_15hpostIR_5do_Rep1 | WT_YR_15hpostIR_5do_Rep2 | WT_YR_12hpostIR_5do_Rep1 | WT_YR_12hpostIR_5do_Rep2 | WT_YR_24hpostIR_6do_Rep1 | WT_YR_24hpostIR_6do_Rep2 | WT_YR_3hpostIR_5do_Rep1 | WT_YR_3hpostIR_5do_Rep2 | WT_YR_6hpostIR_5do_Rep1 | WT_YR_6hpostIR_5do_Rep2 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Sample.LongId | Atha.Ws.sl.ATM1.HZE.1.5h.postIR.5do.Rep1.Array | Atha.Ws.sl.ATM1.HZE.1.5h.postIR.5do.Rep2.Array | Atha.Ws.sl.ATM1.HZE.12h.postIR.5do.Rep1.Array | Atha.Ws.sl.ATM1.HZE.12h.postIR.5do.Rep2.Array | Atha.Ws.sl.ATM1.HZE.24h.postIR.6do.Rep1.Array | Atha.Ws.sl.ATM1.HZE.24h.postIR.6do.Rep2.Array | Atha.Ws.sl.ATM1.HZE.6h.postIR.5do.Rep1.Array | Atha.Ws.sl.ATM1.HZE.6h.postIR.5do.Rep2.Array | Atha.Ws.sl.ATM1.noIR.CTRL.5do.Rep1.Array | Atha.Ws.sl.ATM1.noIR.CTRL.5do.Rep2.Array | Atha.Ws.sl.ATM1.noIR.CTRL.6do.Rep1.Array | Atha.Ws.sl.ATM1.noIR.CTRL.6do.Rep2.Array | Atha.Ws.sl.ATM1.YR.1.5h.postIR.5do.Rep1.Array | Atha.Ws.sl.ATM1.YR.1.5h.postIR.5do.Rep2.Array | Atha.Ws.sl.ATM1.YR.12h.postIR.5do.Rep1.Array | Atha.Ws.sl.ATM1.YR.12h.postIR.5do.Rep2.Array | Atha.Ws.sl.ATM1.YR.24h.postIR.6do.Rep1.Array | Atha.Ws.sl.ATM1.YR.24h.postIR.6do.Rep2.Array | Atha.Ws.sl.ATM1.YR.6h.postIR.5do.Rep1.Array | Atha.Ws.sl.ATM1.YR.6h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.HZE.1.5h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.HZE.1.5h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.HZE.12h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.HZE.12h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.HZE.24h.postIR.6do.Rep1.Array | Atha.Ws.sl.WT.HZE.24h.postIR.6do.Rep2.Array | Atha.Ws.sl.WT.HZE.3h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.HZE.3h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.HZE.6h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.HZE.6h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.noIR.CTRL.5do.Rep1.Array | Atha.Ws.sl.WT.noIR.CTRL.5do.Rep2.Array | Atha.Ws.sl.WT.noIR.CTRL.6do.Rep1.Array | Atha.Ws.sl.WT.noIR.CTRL.6do.Rep2.Array | Atha.Ws.sl.WT.YR.1.5h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.YR.1.5h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.YR.12h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.YR.12h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.YR.24h.postIR.6do.Rep1.Array | Atha.Ws.sl.WT.YR.24h.postIR.6do.Rep2.Array | Atha.Ws.sl.WT.YR.3h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.YR.3h.postIR.5do.Rep2.Array | Atha.Ws.sl.WT.YR.6h.postIR.5do.Rep1.Array | Atha.Ws.sl.WT.YR.6h.postIR.5do.Rep2.Array |
| Sample.Id | Atha.Ws.sl.ATM1.HZE.1.5h.postIR.5do_Rep1 | Atha.Ws.sl.ATM1.HZE.1.5h.postIR.5do_Rep2 | Atha.Ws.sl.ATM1.HZE.12h.postIR.5do_Rep1 | Atha.Ws.sl.ATM1.HZE.12h.postIR.5do_Rep2 | Atha.Ws.sl.ATM1.HZE.24h.postIR.6do_Rep1 | Atha.Ws.sl.ATM1.HZE.24h.postIR.6do_Rep2 | Atha.Ws.sl.ATM1.HZE.6h.postIR.5do_Rep1 | Atha.Ws.sl.ATM1.HZE.6h.postIR.5do_Rep2 | Atha.Ws.sl.ATM1.noIR.CTRL.5do_Rep1 | Atha.Ws.sl.ATM1.noIR.CTRL.5do_Rep2 | Atha.Ws.sl.ATM1.noIR.CTRL.6do_Rep1 | Atha.Ws.sl.ATM1.noIR.CTRL.6do_Rep2 | Atha.Ws.sl.ATM1.YR.1.5h.postIR.5do_Rep1 | Atha.Ws.sl.ATM1.YR.1.5h.postIR.5do_Rep2 | Atha.Ws.sl.ATM1.YR.12h.postIR.5do_Rep1 | Atha.Ws.sl.ATM1.YR.12h.postIR.5do_Rep2 | Atha.Ws.sl.ATM1.YR.24h.postIR.6do_Rep1 | Atha.Ws.sl.ATM1.YR.24h.postIR.6do_Rep2 | Atha.Ws.sl.ATM1.YR.6h.postIR.5do_Rep1 | Atha.Ws.sl.ATM1.YR.6h.postIR.5do_Rep2 | Atha.Ws.sl.WT.HZE.1.5h.postIR.5do_Rep1 | Atha.Ws.sl.WT.HZE.1.5h.postIR.5do_Rep2 | Atha.Ws.sl.WT.HZE.12h.postIR.5do_Rep1 | Atha.Ws.sl.WT.HZE.12h.postIR.5do_Rep2 | Atha.Ws.sl.WT.HZE.24h.postIR.6do_Rep1 | Atha.Ws.sl.WT.HZE.24h.postIR.6do_Rep2 | Atha.Ws.sl.WT.HZE.3h.postIR.5do_Rep1 | Atha.Ws.sl.WT.HZE.3h.postIR.5do_Rep2 | Atha.Ws.sl.WT.HZE.6h.postIR.5do_Rep1 | Atha.Ws.sl.WT.HZE.6h.postIR.5do_Rep2 | Atha.Ws.sl.WT.noIR.CTRL.5do_Rep1 | Atha.Ws.sl.WT.noIR.CTRL.5do_Rep2 | Atha.Ws.sl.WT.noIR.CTRL.6do_Rep1 | Atha.Ws.sl.WT.noIR.CTRL.6do_Rep2 | Atha.Ws.sl.WT.YR.1.5h.postIR.5do_Rep1 | Atha.Ws.sl.WT.YR.1.5h.postIR.5do_Rep2 | Atha.Ws.sl.WT.YR.12h.postIR.5do_Rep1 | Atha.Ws.sl.WT.YR.12h.postIR.5do_Rep2 | Atha.Ws.sl.WT.YR.24h.postIR.6do_Rep1 | Atha.Ws.sl.WT.YR.24h.postIR.6do_Rep2 | Atha.Ws.sl.WT.YR.3h.postIR.5do_Rep1 | Atha.Ws.sl.WT.YR.3h.postIR.5do_Rep2 | Atha.Ws.sl.WT.YR.6h.postIR.5do_Rep1 | Atha.Ws.sl.WT.YR.6h.postIR.5do_Rep2 |
| Sample.Name | Atha_Ws_sl_ATM1_HZE_1.5h-postIR_5do_Rep1 | Atha_Ws_sl_ATM1_HZE_1.5h-postIR_5do_Rep2 | Atha_Ws_sl_ATM1_HZE_12h-postIR_5do_Rep1 | Atha_Ws_sl_ATM1_HZE_12h-postIR_5do_Rep2 | Atha_Ws_sl_ATM1_HZE_24h-postIR_6do_Rep1 | Atha_Ws_sl_ATM1_HZE_24h-postIR_6do_Rep2 | Atha_Ws_sl_ATM1_HZE_6h-postIR_5do_Rep1 | Atha_Ws_sl_ATM1_HZE_6h-postIR_5do_Rep2 | Atha_Ws_sl_ATM1_noIR_CTRL_5do_Rep1 | Atha_Ws_sl_ATM1_noIR_CTRL_5do_Rep2 | Atha_Ws_sl_ATM1_noIR_CTRL_6do_Rep1 | Atha_Ws_sl_ATM1_noIR_CTRL_6do_Rep2 | Atha_Ws_sl_ATM1_YR_1.5h-postIR_5do_Rep1 | Atha_Ws_sl_ATM1_YR_1.5h-postIR_5do_Rep2 | Atha_Ws_sl_ATM1_YR_12h-postIR_5do_Rep1 | Atha_Ws_sl_ATM1_YR_12h-postIR_5do_Rep2 | Atha_Ws_sl_ATM1_YR_24h-postIR_6do_Rep1 | Atha_Ws_sl_ATM1_YR_24h-postIR_6do_Rep2 | Atha_Ws_sl_ATM1_YR_6h-postIR_5do_Rep1 | Atha_Ws_sl_ATM1_YR_6h-postIR_5do_Rep2 | Atha_Ws_sl_WT_HZE_1.5h-postIR_5do_Rep1 | Atha_Ws_sl_WT_HZE_1.5h-postIR_5do_Rep2 | Atha_Ws_sl_WT_HZE_12h-postIR_5do_Rep1 | Atha_Ws_sl_WT_HZE_12h-postIR_5do_Rep2 | Atha_Ws_sl_WT_HZE_24h-postIR_6do_Rep1 | Atha_Ws_sl_WT_HZE_24h-postIR_6do_Rep2 | Atha_Ws_sl_WT_HZE_3h-postIR_5do_Rep1 | Atha_Ws_sl_WT_HZE_3h-postIR_5do_Rep2 | Atha_Ws_sl_WT_HZE_6h-postIR_5do_Rep1 | Atha_Ws_sl_WT_HZE_6h-postIR_5do_Rep2 | Atha_Ws_sl_WT_noIR_CTRL_5do_Rep1 | Atha_Ws_sl_WT_noIR_CTRL_5do_Rep2 | Atha_Ws_sl_WT_noIR_CTRL_6do_Rep1 | Atha_Ws_sl_WT_noIR_CTRL_6do_Rep2 | Atha_Ws_sl_WT_YR_1.5h-postIR_5do_Rep1 | Atha_Ws_sl_WT_YR_1.5h-postIR_5do_Rep2 | Atha_Ws_sl_WT_YR_12h-postIR_5do_Rep1 | Atha_Ws_sl_WT_YR_12h-postIR_5do_Rep2 | Atha_Ws_sl_WT_YR_24h-postIR_6do_Rep1 | Atha_Ws_sl_WT_YR_24h-postIR_6do_Rep2 | Atha_Ws_sl_WT_YR_3h-postIR_5do_Rep1 | Atha_Ws_sl_WT_YR_3h-postIR_5do_Rep2 | Atha_Ws_sl_WT_YR_6h-postIR_5do_Rep1 | Atha_Ws_sl_WT_YR_6h-postIR_5do_Rep2 |
| GLDS | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 | 46 |
| Accession | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 | GLDS-46 |
| Hardware | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish | Petri dish |
| Tissue | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings | Seedlings |
| Age | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days | 5 days |
| Organism | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana | Arabidopsis thaliana |
| Ecotype | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 | WS-0 |
| Genotype | atm1-1 | atm1-1 | atm1-1 | atm1-1 | atm1-1 | atm1-1 | WT | WT | atm1-1 | WT | atm1-1 | atm1-1 | WT | WT | atm1-1 | WT | atm1-1 | atm1-1 | WT | WT | atm1-1 | atm1-1 | WT | atm1-1 | WT | atm1-1 | WT | WT | atm1-1 | WT | WT | WT | atm1-1 | WT | atm1-1 | WT | WT | WT | WT | atm1-1 | WT | WT | WT | WT |
| Variety | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 atm1-1 | WS-0 atm1-1 | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 WT | WS-0 WT | WS-0 atm1-1 | WS-0 WT | WS-0 WT | WS-0 WT | WS-0 WT |
| Radiation | GA6.0. | GA1.5. | HA1.5 | HA12.0 | HA24.0 | GA1.5 | GW12.0 | HW12.0 | GA6.0 | GW6.0 | GA12.0 | GA24.0 | HW12.0 | GW3.0 | HA24.0 | GW12.0 | GA24.0 | HA1.5 | HW6.0 | HW1.5 | UA | UA | GW24.0 | UA.day.2 | HW24.0 | GA12.0 | GW1.5 | GW1.5 | HA6.0 | HW1.5 | UW.day.2 | HW3.0 | HA12.0 | HW3.0 | HA6.0 | HW6.0 | HW24.0 | GW6.0 | UW | UA | UW | GW3.0 | GW24.0 | UW |
| Gravity | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial | Terrestrial |
| Developmental | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. | 5 day old seedling, cotyledons opened fully. |
| Time.series.or.Concentration.gradient | 1.5h-postIR | 1.5h-postIR | 12h-postIR | 12h-postIR | 24h-postIR | 24h-postIR | 6h-postIR | 6h-postIR | NoIR | NoIR | NoIR | NoIR | 1.5h-postIR | 1.5h-postIR | 12h-postIR | 12h-postIR | 24h-postIR | 24h-postIR | 6h-postIR | 6h-postIR | 1.5h-postIR | 1.5h-postIR | 12h-postIR | 12h-postIR | 24h-postIR | 24h-postIR | 3h-postIR | 3h-postIR | 6h-postIR | 6h-postIR | NoIR | NoIR | NoIR | NoIR | 1.5h-postIR | 1.5h-postIR | 12h-postIR | 12h-postIR | 24h-postIR | 24h-postIR | 3h-postIR | 3h-postIR | 6h-postIR | 6h-postIR |
| Light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light | Cool white light |
| Assay..RNAseq. | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling | Microarray Transcription Profiling |
| Temperature | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 | 20-22 |
| Treatment.type | HZE | HZE | HZE | HZE | HZE | HZE | HZE | HZE | control | control | control | control | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | HZE | HZE | HZE | HZE | HZE | HZE | HZE | HZE | HZE | HZE | control | control | control | control | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma | Gamma |
| Treatment.intensity | GA6.0. | GA1.5. | HA1.5 | HA12.0 | HA24.0 | GA1.5 | GW12.0 | HW12.0 | GA6.0 | GW6.0 | GA12.0 | GA24.0 | HW12.0 | GW3.0 | HA24.0 | GW12.0 | GA24.0 | HA1.5 | HW6.0 | HW1.5 | UA | UA | GW24.0 | UA.day.2 | HW24.0 | GA12.0 | GW1.5 | GW1.5 | HA6.0 | HW1.5 | UW.day.2 | HW3.0 | HA12.0 | HW3.0 | HA6.0 | HW6.0 | HW24.0 | GW6.0 | UW | UA | UW | GW3.0 | GW24.0 | UW |
| Treament.timing | seedling_6h after Gamma irradiation_atm-1 | seedling_1.5h after Gamma irradiation_atm-1 | seedling_1.5h after HZE irradiation_atm-1 | seedling_12h after HZE irradiation_atm-1 | seedling_24h after HZE irradiation_atm-1 | seedling_1.5h after Gamma irradiation_atm-1 | seedling_12h after Gamma irradiation_WT | seedling_12h after HZE irradiation_WT | seedling_6h after Gamma irradiation_atm-1 | seedling_6h after Gamma irradiation_WT | seedling_12h after Gamma irradiation_atm-1 | seedling_24h after Gamma irradiation_atm-1 | seedling_12h after HZE irradiation_WT | seedling_3h after Gamma irradiation_WT | seedling_24h after HZE irradiation_atm-1 | seedling_12h after Gamma irradiation_WT | seedling_24h after Gamma irradiation_atm-1 | seedling_1.5h after HZE irradiation_atm-1 | seedling_6h after HZE irradiation_WT | seedling_1.5h after HZE irradiation_WT | seedling_day 1_unirradiated control_atm-1 | seedling_day 1_unirradiated control_atm-1 | seedling_24h after Gamma irradiation_WT | seedling_day 2_unirradiated control_atm-1 | seedling_24h after HZE irradiation_WT | seedling_12h after Gamma irradiation_atm-1 | seedling_1.5h after Gamma irradiation_WT | seedling_1.5h after Gamma irradiation_WT | seedling_6h after HZE irradiation_atm-1 | seedling_1.5h after HZE irradiation_WT | seedling_day 2_unirradiated control_WT | seedling_3h after HZE irradiation_WT | seedling_12h after HZE irradiation_atm-1 | seedling_3h after HZE irradiation_WT | seedling_6h after HZE irradiation_atm-1 | seedling_6h after HZE irradiation_WT | seedling_24h after HZE irradiation_WT | seedling_6h after Gamma irradiation_WT | seedling_day 2_unirradiated control_WT | seedling_day 2_unirradiated control_atm-1 | seedling_day 1_unirradiated control_WT | seedling_3h after Gamma irradiation_WT | seedling_24h after Gamma irradiation_WT | seedling_day 1_unirradiated control_WT |
| Preservation.Method. | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |
readData.out <- readData(inputFile)
## Warning in DESeqDataSet(se, design = design, ignoreRank): some variables in
## design formula are characters, converting to factors
kable( head(readData.out$data) ) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| ATM1_HZE_15hpostIR_5do_Rep1 | ATM1_HZE_15hpostIR_5do_Rep2 | ATM1_HZE_12hpostIR_5do_Rep1 | ATM1_HZE_12hpostIR_5do_Rep2 | ATM1_HZE_24hpostIR_6do_Rep1 | ATM1_HZE_24hpostIR_6do_Rep2 | ATM1_HZE_6hpostIR_5do_Rep1 | ATM1_HZE_6hpostIR_5do_Rep2 | ATM1_noIR_CTRL_5do_Rep1 | ATM1_noIR_CTRL_5do_Rep2 | ATM1_noIR_CTRL_6do_Rep1 | ATM1_noIR_CTRL_6do_Rep2 | ATM1_YR_15hpostIR_5do_Rep1 | ATM1_YR_15hpostIR_5do_Rep2 | ATM1_YR_12hpostIR_5do_Rep1 | ATM1_YR_12hpostIR_5do_Rep2 | ATM1_YR_24hpostIR_6do_Rep1 | ATM1_YR_24hpostIR_6do_Rep2 | ATM1_YR_6hpostIR_5do_Rep1 | ATM1_YR_6hpostIR_5do_Rep2 | WT_HZE_15hpostIR_5do_Rep1 | WT_HZE_15hpostIR_5do_Rep2 | WT_HZE_12hpostIR_5do_Rep1 | WT_HZE_12hpostIR_5do_Rep2 | WT_HZE_24hpostIR_6do_Rep1 | WT_HZE_24hpostIR_6do_Rep2 | WT_HZE_3hpostIR_5do_Rep1 | WT_HZE_3hpostIR_5do_Rep2 | WT_HZE_6hpostIR_5do_Rep1 | WT_HZE_6hpostIR_5do_Rep2 | WT_noIR_CTRL_5do_Rep1 | WT_noIR_CTRL_5do_Rep2 | WT_noIR_CTRL_6do_Rep1 | WT_noIR_CTRL_6do_Rep2 | WT_YR_15hpostIR_5do_Rep1 | WT_YR_15hpostIR_5do_Rep2 | WT_YR_12hpostIR_5do_Rep1 | WT_YR_12hpostIR_5do_Rep2 | WT_YR_24hpostIR_6do_Rep1 | WT_YR_24hpostIR_6do_Rep2 | WT_YR_3hpostIR_5do_Rep1 | WT_YR_3hpostIR_5do_Rep2 | WT_YR_6hpostIR_5do_Rep1 | WT_YR_6hpostIR_5do_Rep2 | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| AT2G46830 | 3.700440 | 3.700440 | 3.000000 | 3.000000 | 3.807355 | 3.807355 | 3.321928 | 3.321928 | 3.807355 | 3.807355 | 3.807355 | 3.806773 | 3.807355 | 3.906891 | 3.000000 | 3.000000 | 3.906891 | 3.906891 | 3.599633 | 3.584963 | 3.700440 | 3.700440 | 3.000000 | 3.000000 | 3.700440 | 3.807355 | 3.584963 | 3.584963 | 3.169925 | 3.169925 | 3.906891 | 3.906891 | 3.807355 | 3.807355 | 3.906891 | 3.906891 | 3.000000 | 3.000000 | 3.906891 | 3.906891 | 3.807355 | 3.807355 | 3.700440 | 3.700440 |
| AT4G21070 | 3.169925 | 3.169925 | 3.584963 | 3.584963 | 3.584963 | 3.459432 | 3.459432 | 3.459432 | 3.169925 | 3.000000 | 3.169925 | 3.169472 | 3.584963 | 3.584963 | 3.584963 | 3.584963 | 3.459432 | 3.459432 | 3.599633 | 3.584963 | 4.000000 | 4.000000 | 3.906891 | 3.906891 | 3.700440 | 3.700440 | 4.000000 | 4.000000 | 3.906891 | 3.906891 | 3.000000 | 3.000000 | 3.169925 | 3.169925 | 4.000000 | 4.000000 | 3.700440 | 3.700440 | 3.584963 | 3.584963 | 3.906891 | 3.906891 | 3.807355 | 3.807355 |
| AT1G32900 | 3.584963 | 3.584963 | 3.169925 | 3.169925 | 3.700440 | 3.807355 | 3.700440 | 3.700440 | 3.906891 | 3.906891 | 3.807355 | 3.806773 | 4.000000 | 4.000000 | 3.000000 | 3.000000 | 4.000000 | 4.000000 | 3.599633 | 3.700440 | 3.584963 | 3.584963 | 3.169925 | 3.321928 | 3.700440 | 3.700440 | 3.459432 | 3.459432 | 3.584963 | 3.584963 | 3.906891 | 3.906891 | 3.700440 | 3.700440 | 4.000000 | 4.000000 | 3.169925 | 3.000000 | 4.000000 | 4.000000 | 3.906891 | 3.906891 | 3.700440 | 3.584963 |
| AT3G02380 | 3.584963 | 3.584963 | 3.169925 | 3.169925 | 3.807355 | 3.807355 | 3.459432 | 3.321928 | 3.807355 | 3.807355 | 3.807355 | 3.806773 | 3.906891 | 3.906891 | 3.169925 | 3.169925 | 3.906891 | 3.906891 | 3.715671 | 3.700440 | 3.584963 | 3.584963 | 3.169925 | 3.169925 | 3.807355 | 3.807355 | 3.321928 | 3.321928 | 3.169925 | 3.169925 | 3.906891 | 3.906891 | 3.807355 | 3.807355 | 3.906891 | 3.906891 | 3.169925 | 3.169925 | 3.807355 | 3.906891 | 3.807355 | 3.807355 | 3.700440 | 3.700440 |
| AT1G10070 | 4.000000 | 4.000000 | 3.169925 | 3.321928 | 3.321928 | 3.459432 | 3.459432 | 3.459432 | 3.321928 | 3.321928 | 3.459432 | 3.458913 | 3.169925 | 3.169925 | 3.169925 | 3.459432 | 3.169925 | 3.169925 | 3.182161 | 3.321928 | 4.000000 | 4.000000 | 3.169925 | 3.169925 | 3.169925 | 3.321928 | 4.000000 | 4.000000 | 3.321928 | 3.459432 | 3.169925 | 3.169925 | 3.459432 | 3.321928 | 3.169925 | 3.169925 | 3.169925 | 3.169925 | 3.169925 | 3.169925 | 3.169925 | 3.321928 | 3.169925 | 3.459432 |
| AT5G48720 | 3.459432 | 3.459432 | 3.700440 | 3.700440 | 3.584963 | 3.584963 | 3.584963 | 3.584963 | 3.321928 | 3.321928 | 3.321928 | 3.321439 | 3.584963 | 3.700440 | 3.700440 | 3.700440 | 3.584963 | 3.584963 | 3.599633 | 3.700440 | 4.087463 | 4.087463 | 4.000000 | 3.906891 | 3.807355 | 3.807355 | 4.000000 | 4.000000 | 4.000000 | 4.000000 | 3.321928 | 3.321928 | 3.321928 | 3.321928 | 4.087463 | 4.087463 | 3.807355 | 3.807355 | 3.700440 | 3.700440 | 4.000000 | 4.000000 | 3.906891 | 3.906891 |
readSampleInfo.out <- readSampleInfo(sampleInfoFile)
kable( readSampleInfo.out ) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| Genotype | Radiation | Time.series.or.Concentration.gradient | Treatment.type | Treatment.intensity | Treament.timing | |
|---|---|---|---|---|---|---|
| ATM1_HZE_15hpostIR_5do_Rep1 | Genotypeatm11 | RadiationGA60 | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeHZE | Treatment.intensityGA60 | Treament.timingseedling_6h after Gamma irradiation_atm1 |
| ATM1_HZE_15hpostIR_5do_Rep2 | Genotypeatm11 | RadiationGA15 | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeHZE | Treatment.intensityGA15 | Treament.timingseedling_15h after Gamma irradiation_atm1 |
| ATM1_HZE_12hpostIR_5do_Rep1 | Genotypeatm11 | RadiationHA15 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeHZE | Treatment.intensityHA15 | Treament.timingseedling_15h after HZE irradiation_atm1 |
| ATM1_HZE_12hpostIR_5do_Rep2 | Genotypeatm11 | RadiationHA120 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeHZE | Treatment.intensityHA120 | Treament.timingseedling_12h after HZE irradiation_atm1 |
| ATM1_HZE_24hpostIR_6do_Rep1 | Genotypeatm11 | RadiationHA240 | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeHZE | Treatment.intensityHA240 | Treament.timingseedling_24h after HZE irradiation_atm1 |
| ATM1_HZE_24hpostIR_6do_Rep2 | Genotypeatm11 | RadiationGA15 | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeHZE | Treatment.intensityGA15 | Treament.timingseedling_15h after Gamma irradiation_atm1 |
| ATM1_HZE_6hpostIR_5do_Rep1 | GenotypeWT | RadiationGW120 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeHZE | Treatment.intensityGW120 | Treament.timingseedling_12h after Gamma irradiation_WT |
| ATM1_HZE_6hpostIR_5do_Rep2 | GenotypeWT | RadiationHW120 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeHZE | Treatment.intensityHW120 | Treament.timingseedling_12h after HZE irradiation_WT |
| ATM1_noIR_CTRL_5do_Rep1 | Genotypeatm11 | RadiationGA60 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityGA60 | Treament.timingseedling_6h after Gamma irradiation_atm1 |
| ATM1_noIR_CTRL_5do_Rep2 | GenotypeWT | RadiationGW60 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityGW60 | Treament.timingseedling_6h after Gamma irradiation_WT |
| ATM1_noIR_CTRL_6do_Rep1 | Genotypeatm11 | RadiationGA120 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityGA120 | Treament.timingseedling_12h after Gamma irradiation_atm1 |
| ATM1_noIR_CTRL_6do_Rep2 | Genotypeatm11 | RadiationGA240 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityGA240 | Treament.timingseedling_24h after Gamma irradiation_atm1 |
| ATM1_YR_15hpostIR_5do_Rep1 | GenotypeWT | RadiationHW120 | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeGamma | Treatment.intensityHW120 | Treament.timingseedling_12h after HZE irradiation_WT |
| ATM1_YR_15hpostIR_5do_Rep2 | GenotypeWT | RadiationGW30 | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeGamma | Treatment.intensityGW30 | Treament.timingseedling_3h after Gamma irradiation_WT |
| ATM1_YR_12hpostIR_5do_Rep1 | Genotypeatm11 | RadiationHA240 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeGamma | Treatment.intensityHA240 | Treament.timingseedling_24h after HZE irradiation_atm1 |
| ATM1_YR_12hpostIR_5do_Rep2 | GenotypeWT | RadiationGW120 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeGamma | Treatment.intensityGW120 | Treament.timingseedling_12h after Gamma irradiation_WT |
| ATM1_YR_24hpostIR_6do_Rep1 | Genotypeatm11 | RadiationGA240 | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeGamma | Treatment.intensityGA240 | Treament.timingseedling_24h after Gamma irradiation_atm1 |
| ATM1_YR_24hpostIR_6do_Rep2 | Genotypeatm11 | RadiationHA15 | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeGamma | Treatment.intensityHA15 | Treament.timingseedling_15h after HZE irradiation_atm1 |
| ATM1_YR_6hpostIR_5do_Rep1 | GenotypeWT | RadiationHW60 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeGamma | Treatment.intensityHW60 | Treament.timingseedling_6h after HZE irradiation_WT |
| ATM1_YR_6hpostIR_5do_Rep2 | GenotypeWT | RadiationHW15 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeGamma | Treatment.intensityHW15 | Treament.timingseedling_15h after HZE irradiation_WT |
| WT_HZE_15hpostIR_5do_Rep1 | Genotypeatm11 | RadiationUA | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeHZE | Treatment.intensityUA | Treament.timingseedling_day 1_unirradiated control_atm1 |
| WT_HZE_15hpostIR_5do_Rep2 | Genotypeatm11 | RadiationUA | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeHZE | Treatment.intensityUA | Treament.timingseedling_day 1_unirradiated control_atm1 |
| WT_HZE_12hpostIR_5do_Rep1 | GenotypeWT | RadiationGW240 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeHZE | Treatment.intensityGW240 | Treament.timingseedling_24h after Gamma irradiation_WT |
| WT_HZE_12hpostIR_5do_Rep2 | Genotypeatm11 | RadiationUAday2 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeHZE | Treatment.intensityUAday2 | Treament.timingseedling_day 2_unirradiated control_atm1 |
| WT_HZE_24hpostIR_6do_Rep1 | GenotypeWT | RadiationHW240 | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeHZE | Treatment.intensityHW240 | Treament.timingseedling_24h after HZE irradiation_WT |
| WT_HZE_24hpostIR_6do_Rep2 | Genotypeatm11 | RadiationGA120 | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeHZE | Treatment.intensityGA120 | Treament.timingseedling_12h after Gamma irradiation_atm1 |
| WT_HZE_3hpostIR_5do_Rep1 | GenotypeWT | RadiationGW15 | Time.series.or.Concentration.gradient3hpostIR | Treatment.typeHZE | Treatment.intensityGW15 | Treament.timingseedling_15h after Gamma irradiation_WT |
| WT_HZE_3hpostIR_5do_Rep2 | GenotypeWT | RadiationGW15 | Time.series.or.Concentration.gradient3hpostIR | Treatment.typeHZE | Treatment.intensityGW15 | Treament.timingseedling_15h after Gamma irradiation_WT |
| WT_HZE_6hpostIR_5do_Rep1 | Genotypeatm11 | RadiationHA60 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeHZE | Treatment.intensityHA60 | Treament.timingseedling_6h after HZE irradiation_atm1 |
| WT_HZE_6hpostIR_5do_Rep2 | GenotypeWT | RadiationHW15 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeHZE | Treatment.intensityHW15 | Treament.timingseedling_15h after HZE irradiation_WT |
| WT_noIR_CTRL_5do_Rep1 | GenotypeWT | RadiationUWday2 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityUWday2 | Treament.timingseedling_day 2_unirradiated control_WT |
| WT_noIR_CTRL_5do_Rep2 | GenotypeWT | RadiationHW30 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityHW30 | Treament.timingseedling_3h after HZE irradiation_WT |
| WT_noIR_CTRL_6do_Rep1 | Genotypeatm11 | RadiationHA120 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityHA120 | Treament.timingseedling_12h after HZE irradiation_atm1 |
| WT_noIR_CTRL_6do_Rep2 | GenotypeWT | RadiationHW30 | Time.series.or.Concentration.gradientNoIR | Treatment.typecontrol | Treatment.intensityHW30 | Treament.timingseedling_3h after HZE irradiation_WT |
| WT_YR_15hpostIR_5do_Rep1 | Genotypeatm11 | RadiationHA60 | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeGamma | Treatment.intensityHA60 | Treament.timingseedling_6h after HZE irradiation_atm1 |
| WT_YR_15hpostIR_5do_Rep2 | GenotypeWT | RadiationHW60 | Time.series.or.Concentration.gradient15hpostIR | Treatment.typeGamma | Treatment.intensityHW60 | Treament.timingseedling_6h after HZE irradiation_WT |
| WT_YR_12hpostIR_5do_Rep1 | GenotypeWT | RadiationHW240 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeGamma | Treatment.intensityHW240 | Treament.timingseedling_24h after HZE irradiation_WT |
| WT_YR_12hpostIR_5do_Rep2 | GenotypeWT | RadiationGW60 | Time.series.or.Concentration.gradient12hpostIR | Treatment.typeGamma | Treatment.intensityGW60 | Treament.timingseedling_6h after Gamma irradiation_WT |
| WT_YR_24hpostIR_6do_Rep1 | GenotypeWT | RadiationUW | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeGamma | Treatment.intensityUW | Treament.timingseedling_day 2_unirradiated control_WT |
| WT_YR_24hpostIR_6do_Rep2 | Genotypeatm11 | RadiationUA | Time.series.or.Concentration.gradient24hpostIR | Treatment.typeGamma | Treatment.intensityUA | Treament.timingseedling_day 2_unirradiated control_atm1 |
| WT_YR_3hpostIR_5do_Rep1 | GenotypeWT | RadiationUW | Time.series.or.Concentration.gradient3hpostIR | Treatment.typeGamma | Treatment.intensityUW | Treament.timingseedling_day 1_unirradiated control_WT |
| WT_YR_3hpostIR_5do_Rep2 | GenotypeWT | RadiationGW30 | Time.series.or.Concentration.gradient3hpostIR | Treatment.typeGamma | Treatment.intensityGW30 | Treament.timingseedling_3h after Gamma irradiation_WT |
| WT_YR_6hpostIR_5do_Rep1 | GenotypeWT | RadiationGW240 | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeGamma | Treatment.intensityGW240 | Treament.timingseedling_24h after Gamma irradiation_WT |
| WT_YR_6hpostIR_5do_Rep2 | GenotypeWT | RadiationUW | Time.series.or.Concentration.gradient6hpostIR | Treatment.typeGamma | Treatment.intensityUW | Treament.timingseedling_day 1_unirradiated control_WT |
input_selectOrg ="NEW"
input_selectGO <- 'GOBP' #Gene set category
input_noIDConversion = TRUE
allGeneInfo.out <- geneInfo(geneInfoFile)
converted.out = NULL
convertedData.out <- convertedData()
nGenesFilter()
## [1] "16156 genes in 44 samples. 16156 genes passed filter.\n Original gene IDs used."
convertedCounts.out <- convertedCounts() # converted counts, just for compatibility
# Read counts per library
parDefault = par()
par(mar=c(12,4,2,2))
# barplot of total read counts
x <- readData.out$rawCounts
groups = as.factor( detectGroups(colnames(x ) ) )
if(nlevels(groups)<=1 | nlevels(groups) >20 )
col1 = 'green' else
col1 = rainbow(nlevels(groups))[ groups ]
barplot( colSums(x)/1e6,
col=col1,las=3, main="Total read counts (millions)")
readCountsBias() # detecting bias in sequencing depth
## [1] 0.1982081
## [1] 0.9367154
## [1] 0.5100945
## [1] 0.6402451
## [1] 0.2308868
## [1] 0.5100945
## [1] 0.527478
## [1] "No bias detected"
# Box plot
x = readData.out$data
boxplot(x, las = 2, col=col1,
ylab='Transformed expression levels',
main='Distribution of transformed data')
#Density plot
par(parDefault)
## Warning in par(parDefault): graphical parameter "cin" cannot be set
## Warning in par(parDefault): graphical parameter "cra" cannot be set
## Warning in par(parDefault): graphical parameter "csi" cannot be set
## Warning in par(parDefault): graphical parameter "cxy" cannot be set
## Warning in par(parDefault): graphical parameter "din" cannot be set
## Warning in par(parDefault): graphical parameter "page" cannot be set
densityPlot()
# Scatter plot of the first two samples
plot(x[,1:2],xlab=colnames(x)[1],ylab=colnames(x)[2],
main='Scatter plot of first two samples')
####plot gene or gene family
input_selectOrg ="BestMatch"
input_geneSearch <- 'HOXA' #Gene ID for searching
genePlot()
## NULL
input_useSD <- 'FALSE' #Use standard deviation instead of standard error in error bar?
geneBarPlotError()
## NULL
# hierarchical clustering tree
x <- readData.out$data
maxGene <- apply(x,1,max)
# remove bottom 25% lowly expressed genes, which inflate the PPC
x <- x[which(maxGene > quantile(maxGene)[1] ) ,]
plot(as.dendrogram(hclust2( dist2(t(x)))), ylab="1 - Pearson C.C.", type = "rectangle")
#Correlation matrix
input_labelPCC <- TRUE #Show correlation coefficient?
correlationMatrix()
# Parameters for heatmap
input_nGenes <- 1000 #Top genes for heatmap
input_geneCentering <- TRUE #centering genes ?
input_sampleCentering <- FALSE #Center by sample?
input_geneNormalize <- FALSE #Normalize by gene?
input_sampleNormalize <- FALSE #Normalize by sample?
input_noSampleClustering <- FALSE #Use original sample order
input_heatmapCutoff <- 4 #Remove outliers beyond number of SDs
input_distFunctions <- 1 #which distant funciton to use
input_hclustFunctions <- 1 #Linkage type
input_heatColors1 <- 1 #Colors
input_selectFactorsHeatmap <- 'Treatment.type' #Sample coloring factors
png('heatmap.png', width = 10, height = 15, units = 'in', res = 300)
staticHeatmap()
dev.off()
## png
## 2
[heatmap] (heatmap.png)
heatmapPlotly() # interactive heatmap using Plotly
input_nGenesKNN <- 2000 #Number of genes fro k-Means
input_nClusters <- 4 #Number of clusters
maxGeneClustering = 12000
input_kmeansNormalization <- 'geneMean' #Normalization
input_KmeansReRun <- 0 #Random seed
distributionSD() #Distribution of standard deviations
KmeansNclusters() #Number of clusters
Kmeans.out = Kmeans() #Running K-means
KmeansHeatmap() #Heatmap for k-Means
#Read gene sets for enrichment analysis
sqlite <- dbDriver('SQLite')
input_selectGO3 <- 'GOBP' #Gene set category
input_minSetSize <- 15 #Min gene set size
input_maxSetSize <- 2000 #Max gene set size
GeneSets.out <-readGeneSets( geneSetFile,
convertedData.out, input_selectGO3,input_selectOrg,
c(input_minSetSize, input_maxSetSize) )
# Alternatively, users can use their own GMT files by
#GeneSets.out <- readGMTRobust('somefile.GMT')
results <- KmeansGO() #Enrichment analysis for k-Means clusters
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| Cluster | adj.Pval | Genes | Pathways |
|---|---|---|---|
| A | 1.35e-29 | 81 | Response to abiotic stimulus |
| 3.16e-23 | 44 | Response to light stimulus | |
| 1.07e-22 | 44 | Response to radiation | |
| 3.39e-14 | 21 | Response to red or far red light | |
| 4.95e-12 | 47 | Response to oxygen-containing compound | |
| 1.12e-11 | 45 | Cellular response to chemical stimulus | |
| 5.04e-11 | 46 | Response to hormone | |
| 8.06e-11 | 46 | Response to endogenous stimulus | |
| 1.44e-10 | 50 | Response to organic substance | |
| 2.63e-10 | 33 | Response to inorganic substance | |
| B | 7.89e-29 | 53 | Secondary metabolic process |
| 2.07e-25 | 118 | Small molecule metabolic process | |
| 6.47e-25 | 129 | Response to abiotic stimulus | |
| 5.15e-24 | 34 | Secondary metabolite biosynthetic process | |
| 3.74e-21 | 86 | Oxoacid metabolic process | |
| 5.70e-21 | 86 | Organic acid metabolic process | |
| 4.84e-18 | 28 | Response to karrikin | |
| 8.22e-18 | 44 | Sulfur compound metabolic process | |
| 8.22e-18 | 31 | Sulfur compound biosynthetic process | |
| 8.26e-17 | 89 | Oxidation-reduction process | |
| C | 1.59e-25 | 106 | Response to abiotic stimulus |
| 5.67e-23 | 84 | Cellular response to chemical stimulus | |
| 3.03e-18 | 89 | Response to organic substance | |
| 2.37e-17 | 29 | Cellular response to decreased oxygen levels | |
| 2.37e-17 | 29 | Cellular response to oxygen levels | |
| 1.68e-16 | 28 | Cellular response to hypoxia | |
| 1.68e-16 | 75 | Response to oxygen-containing compound | |
| 3.77e-16 | 75 | Catabolic process | |
| 4.95e-16 | 29 | Response to decreased oxygen levels | |
| 5.13e-16 | 29 | Response to oxygen levels | |
| D | 2.43e-56 | 136 | Response to abiotic stimulus |
| 5.50e-50 | 44 | Response to chitin | |
| 2.01e-44 | 106 | Response to oxygen-containing compound | |
| 2.06e-44 | 48 | Cellular response to decreased oxygen levels | |
| 2.06e-44 | 48 | Cellular response to oxygen levels | |
| 2.06e-44 | 48 | Cellular response to hypoxia | |
| 1.14e-42 | 101 | Cellular response to chemical stimulus | |
| 6.03e-42 | 48 | Response to hypoxia | |
| 1.18e-41 | 48 | Response to decreased oxygen levels | |
| 1.25e-41 | 90 | Cellular response to stress |
input_seedTSNE <- 0 #Random seed for t-SNE
input_colorGenes <- TRUE #Color genes in t-SNE plot?
tSNEgenePlot() #Plot genes using t-SNE
input_selectFactors <- 'Radiation' #Factor coded by color
input_selectFactors2 <- 'Genotype' #Factor coded by shape
input_tsneSeed2 <- 0 #Random seed for t-SNE
#PCA, MDS and t-SNE plots
PCAplot()
MDSplot()
tSNEplot()
#Read gene sets for pathway analysis using PGSEA on principal components
input_selectGO6 <- 'GOBP'
GeneSets.out <-readGeneSets( geneSetFile,
convertedData.out, input_selectGO6,input_selectOrg,
c(input_minSetSize, input_maxSetSize) )
PCApathway() # Run PGSEA analysis
## Warning: Package 'KEGG.db' is deprecated and will be removed from Bioconductor
## version 3.12
cat( PCA2factor() ) #The correlation between PCs with factors
##
## Correlation between Principal Components (PCs) with factors
## PC1 is correlated with Treatment.type (p=6.94e-04).
## PC2 is correlated with Time.series.or.Concentration.gradient (p=3.65e-13).
## PC3 is correlated with Time.series.or.Concentration.gradient (p=3.49e-04).
input_CountsDEGMethod <- 2 #DESeq2= 3,limma-voom=2,limma-trend=1
input_limmaPval <- 0.1 #FDR cutoff
input_limmaFC <- 2 #Fold-change cutoff
input_selectModelComprions <- c('Treatment.type: Treatment.typecontrol vs. Treatment.typeGamma','Treatment.type: Treatment.typecontrol vs. Treatment.typeHZE') #Selected comparisons
input_selectFactorsModel <- 'Treatment.type' #Selected comparisons
input_selectInteractions <- NULL #Selected comparisons
input_selectBlockFactorsModel <- NULL #Selected comparisons
factorReferenceLevels.out <- c('Treatment.type:Treatment.typeHZE')
limma.out <- limma()
DEG.data.out <- DEG.data()
limma.out$comparisons
## [1] "Treatment.typecontrol-Treatment.typeGamma"
## [2] "Treatment.typecontrol-Treatment.typeHZE"
input_selectComparisonsVenn = limma.out$comparisons[1:3] # use first three comparisons
input_UpDownRegulated <- FALSE #Split up and down regulated genes
vennPlot() # Venn diagram
sigGeneStats() # number of DEGs as figure
sigGeneStatsTable() # number of DEGs as table
## Comparisons
## Treatment.typecontrol-Treatment.typeGamma Treatment.typecontrol-Treatment.typeGamma
## Treatment.typecontrol-Treatment.typeHZE Treatment.typecontrol-Treatment.typeHZE
## Up Down
## Treatment.typecontrol-Treatment.typeGamma 0 0
## Treatment.typecontrol-Treatment.typeHZE 0 0
input_selectContrast = limma.out$comparisons[1] # use first comparisons
selectedHeatmap.data.out <- selectedHeatmap.data()
selectedHeatmap() # heatmap for DEGs in selected comparison
## Error in array(x, c(length(x), 1L), if (!is.null(names(x))) list(names(x), : 'data' must be of a vector type, was 'NULL'
# Save gene lists and data into files
write.csv( selectedHeatmap.data()$genes, 'heatmap.data.csv')
write.csv(DEG.data(),'DEG.data.csv' )
write(AllGeneListsGMT() ,'AllGeneListsGMT.gmt')
input_selectGO2 <- 'GOBP' #Gene set category
geneListData.out <- geneListData()
volcanoPlot()
scatterPlot()
## Error in findContrastSamples(input_selectContrast, colnames(convertedData.out), : object 'c.out' not found
MAplot()
## Error in findContrastSamples(input_selectContrast, colnames(convertedData.out), : object 'c.out' not found
geneListGOTable.out <- geneListGOTable()
# Read pathway data again
GeneSets.out <-readGeneSets( geneSetFile,
convertedData.out, input_selectGO2,input_selectOrg,
c(input_minSetSize, input_maxSetSize) )
input_removeRedudantSets <- TRUE #Remove highly redundant gene sets?
results <- geneListGO() #Enrichment analysis
## Error in if (dim(results1)[2] == 1) return(results1) else {: argument is of length zero
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| Cluster | adj.Pval | Genes | Pathways |
|---|---|---|---|
| A | 1.35e-29 | 81 | Response to abiotic stimulus |
| 3.16e-23 | 44 | Response to light stimulus | |
| 1.07e-22 | 44 | Response to radiation | |
| 3.39e-14 | 21 | Response to red or far red light | |
| 4.95e-12 | 47 | Response to oxygen-containing compound | |
| 1.12e-11 | 45 | Cellular response to chemical stimulus | |
| 5.04e-11 | 46 | Response to hormone | |
| 8.06e-11 | 46 | Response to endogenous stimulus | |
| 1.44e-10 | 50 | Response to organic substance | |
| 2.63e-10 | 33 | Response to inorganic substance | |
| B | 7.89e-29 | 53 | Secondary metabolic process |
| 2.07e-25 | 118 | Small molecule metabolic process | |
| 6.47e-25 | 129 | Response to abiotic stimulus | |
| 5.15e-24 | 34 | Secondary metabolite biosynthetic process | |
| 3.74e-21 | 86 | Oxoacid metabolic process | |
| 5.70e-21 | 86 | Organic acid metabolic process | |
| 4.84e-18 | 28 | Response to karrikin | |
| 8.22e-18 | 44 | Sulfur compound metabolic process | |
| 8.22e-18 | 31 | Sulfur compound biosynthetic process | |
| 8.26e-17 | 89 | Oxidation-reduction process | |
| C | 1.59e-25 | 106 | Response to abiotic stimulus |
| 5.67e-23 | 84 | Cellular response to chemical stimulus | |
| 3.03e-18 | 89 | Response to organic substance | |
| 2.37e-17 | 29 | Cellular response to decreased oxygen levels | |
| 2.37e-17 | 29 | Cellular response to oxygen levels | |
| 1.68e-16 | 28 | Cellular response to hypoxia | |
| 1.68e-16 | 75 | Response to oxygen-containing compound | |
| 3.77e-16 | 75 | Catabolic process | |
| 4.95e-16 | 29 | Response to decreased oxygen levels | |
| 5.13e-16 | 29 | Response to oxygen levels | |
| D | 2.43e-56 | 136 | Response to abiotic stimulus |
| 5.50e-50 | 44 | Response to chitin | |
| 2.01e-44 | 106 | Response to oxygen-containing compound | |
| 2.06e-44 | 48 | Cellular response to decreased oxygen levels | |
| 2.06e-44 | 48 | Cellular response to oxygen levels | |
| 2.06e-44 | 48 | Cellular response to hypoxia | |
| 1.14e-42 | 101 | Cellular response to chemical stimulus | |
| 6.03e-42 | 48 | Response to hypoxia | |
| 1.18e-41 | 48 | Response to decreased oxygen levels | |
| 1.25e-41 | 90 | Cellular response to stress |
STRING-db API access. We need to find the taxonomy id of your species, this used by STRING. First we try to guess the ID based on iDEP’s database. Users can also skip this step and assign NCBI taxonomy id directly by findTaxonomyID.out = 10090 # mouse 10090, human 9606 etc.
STRING10_species = read.csv(STRING10_speciesFile)
ix = grep('Arabidopsis thaliana', STRING10_species$official_name )
findTaxonomyID.out <- STRING10_species[ix,1] # find taxonomyID
findTaxonomyID.out
## [1] 3702
Enrichment analysis using STRING
STRINGdb_geneList.out <- STRINGdb_geneList() #convert gene lists
## Error in names(x) <- value: 'names' attribute [2] must be the same length as the vector [1]
input_STRINGdbGO <- 'Process' #'Process', 'Component', 'Function', 'KEGG', 'Pfam', 'InterPro'
results <- stringDB_GO_enrichmentData() # enrichment using STRING
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
|
PPI network retrieval and analysis
input_nGenesPPI <- 100 #Number of top genes for PPI retrieval and analysis
stringDB_network1(1) #Show PPI network
## Error in stringDB_network1(1): object 'STRINGdb_geneList.out' not found
Generating interactive PPI
write(stringDB_network_link(), 'PPI_results.html') # write results to html file
## Error in stringDB_network_link(): object 'STRINGdb_geneList.out' not found
browseURL('PPI_results.html') # open in browser
input_selectContrast1 = limma.out$comparisons[1]
#input_selectContrast1 = limma.out$comparisons[3] # manually set
input_selectGO <- 'GOBP' #Gene set category
#input_selectGO='custom' # if custom gmt file
input_minSetSize <- 15 #Min size for gene set
input_maxSetSize <- 2000 #Max size for gene set
# Read pathway data again
GeneSets.out <-readGeneSets( geneSetFile,
convertedData.out, input_selectGO,input_selectOrg,
c(input_minSetSize, input_maxSetSize) )
input_pathwayPvalCutoff <- 0.2 #FDR cutoff
input_nPathwayShow <- 30 #Top pathways to show
input_absoluteFold <- FALSE #Use absolute values of fold-change?
input_GenePvalCutoff <- 1 #FDR to remove genes
input_pathwayMethod = 1 # 1 GAGE
gagePathwayData.out <- gagePathwayData() # pathway analysis using GAGE
results <- gagePathwayData.out #Enrichment analysis for k-Means clusters
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| Direction | GAGE analysis: Treatment.typecontrol vs Treatment.typeGamma | statistic | Genes | adj.Pval |
|---|---|---|---|---|
| Down | Ribosome biogenesis | -5.4049 | 343 | 8.9e-05 |
| Ribonucleoprotein complex biogenesis | -5.141 | 438 | 1.6e-04 | |
| NcRNA metabolic process | -4.7065 | 425 | 9.4e-04 | |
| RRNA metabolic process | -4.599 | 244 | 1.3e-03 | |
| RRNA processing | -4.4451 | 239 | 2.1e-03 | |
| NcRNA processing | -4.2518 | 357 | 3.5e-03 | |
| DNA repair | -4.2439 | 314 | 3.5e-03 | |
| Cellular response to DNA damage stimulus | -4.1582 | 337 | 4.4e-03 | |
| DNA replication | -3.8322 | 162 | 1.6e-02 | |
| Double-strand break repair | -3.726 | 118 | 2.5e-02 | |
| RNA modification | -3.5134 | 321 | 4.1e-02 | |
| DNA-dependent DNA replication | -3.4959 | 133 | 4.4e-02 | |
| Energy derivation by oxidation of organic compounds | -3.3904 | 130 | 5.8e-02 | |
| DNA unwinding involved in DNA replication | -3.3888 | 18 | 8.2e-02 | |
| Glycogen metabolic process | -3.3484 | 22 | 7.4e-02 | |
| Energy reserve metabolic process | -3.3484 | 22 | 7.4e-02 | |
| Protein targeting to mitochondrion | -3.3157 | 57 | 6.7e-02 | |
| Protein localization to mitochondrion | -3.3124 | 62 | 6.7e-02 | |
| Establishment of protein localization to mitochondrion | -3.3124 | 62 | 6.7e-02 | |
| Sulfur compound metabolic process | -3.3022 | 319 | 6.7e-02 | |
| Double-strand break repair via homologous recombination | -3.2904 | 78 | 6.7e-02 | |
| Response to toxic substance | -3.2335 | 275 | 6.7e-02 | |
| DNA recombination | -3.17 | 151 | 7.4e-02 | |
| Protein import | -3.1549 | 131 | 7.4e-02 | |
| Recombinational repair | -3.0691 | 83 | 1.0e-01 | |
| DNA geometric change | -3.0136 | 59 | 1.0e-01 | |
| DNA duplex unwinding | -3.0136 | 59 | 1.0e-01 | |
| Protein transmembrane transport | -2.9981 | 112 | 1.0e-01 | |
| Sulfur compound biosynthetic process | -2.9951 | 146 | 1.0e-01 | |
| Mitochondrion organization | -2.975 | 151 | 1.0e-01 |
pathwayListData.out = pathwayListData()
enrichmentPlot(pathwayListData.out, 25 )
enrichmentNetwork(pathwayListData.out )
enrichmentNetworkPlotly(pathwayListData.out)
## Warning: `arrange_()` is deprecated as of dplyr 0.7.0.
## Please use `arrange()` instead.
## See vignette('programming') for more help
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.
input_pathwayMethod = 3 # 1 fgsea
fgseaPathwayData.out <- fgseaPathwayData() #Pathway analysis using fgsea
## Warning in fgsea(pathways = gmt, stats = fold, minSize = input_minSetSize, :
## You are trying to run fgseaSimple. It is recommended to use fgseaMultilevel. To
## run fgseaMultilevel, you need to remove the nperm argument in the fgsea function
## call.
## Warning in preparePathwaysAndStats(pathways, stats, minSize, maxSize, gseaParam, : There are ties in the preranked stats (21.57% of the list).
## The order of those tied genes will be arbitrary, which may produce unexpected results.
results <- fgseaPathwayData.out #Enrichment analysis for k-Means clusters
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| Direction | GSEA analysis: Treatment.typecontrol vs Treatment.typeGamma | NES | Genes | adj.Pval |
|---|---|---|---|---|
| Down | Response to ionizing radiation | -2.4119 | 27 | 7.9e-03 |
| Double-strand break repair | -2.3506 | 118 | 7.9e-03 | |
| Double-strand break repair via homologous recombination | -2.2255 | 78 | 7.9e-03 | |
| Recombinational repair | -2.1395 | 83 | 7.9e-03 | |
| Starch catabolic process | -2.0949 | 15 | 7.9e-03 | |
| DNA repair | -2.0816 | 314 | 7.9e-03 | |
| Reciprocal meiotic recombination | -2.0733 | 46 | 7.9e-03 | |
| Homologous recombination | -2.0733 | 46 | 7.9e-03 | |
| DNA unwinding involved in DNA replication | -2.0675 | 18 | 1.2e-02 | |
| DNA recombination | -2.038 | 151 | 7.9e-03 | |
| Cellular response to DNA damage stimulus | -2.0326 | 337 | 7.9e-03 | |
| DNA replication | -2.0307 | 162 | 7.9e-03 | |
| DNA-dependent DNA replication | -2.0208 | 133 | 7.9e-03 | |
| DNA geometric change | -2.0193 | 59 | 7.9e-03 | |
| DNA duplex unwinding | -2.0193 | 59 | 7.9e-03 | |
| Ribosome biogenesis | -1.9984 | 343 | 7.9e-03 | |
| DNA biosynthetic process | -1.9908 | 57 | 7.9e-03 | |
| Meiosis I | -1.9892 | 56 | 7.9e-03 | |
| Protein targeting to mitochondrion | -1.9849 | 57 | 7.9e-03 | |
| Homologous chromosome segregation | -1.9839 | 26 | 7.9e-03 | |
| Synapsis | -1.9771 | 20 | 1.2e-02 | |
| Glycogen metabolic process | -1.9602 | 22 | 1.2e-02 | |
| Energy reserve metabolic process | -1.9602 | 22 | 1.2e-02 | |
| Toxin catabolic process | -1.9598 | 35 | 1.4e-02 | |
| Protein localization to mitochondrion | -1.9522 | 62 | 1.4e-02 | |
| Establishment of protein localization to mitochondrion | -1.9522 | 62 | 1.4e-02 | |
| Meiosis I cell cycle process | -1.9494 | 62 | 1.4e-02 | |
| Regulation of sulfur metabolic process | -1.9305 | 28 | 1.2e-02 | |
| S-glycoside biosynthetic process | -1.9286 | 37 | 1.7e-02 | |
| Up | Branched-chain amino acid catabolic process | 2.0016 | 16 | 7.9e-03 |
pathwayListData.out = pathwayListData()
enrichmentPlot(pathwayListData.out, 25 )
enrichmentNetwork(pathwayListData.out )
enrichmentNetworkPlotly(pathwayListData.out)
PGSEAplot() # pathway analysis using PGSEA
## Error in findContrastSamples(input_selectContrast1, colnames(convertedData.out), : object 'c.out' not found
input_selectContrast2 = limma.out$comparisons[1]
#input_selectContrast2 = limma.out$comparisons[3] # manually set
input_limmaPvalViz <- 0.1 #FDR to filter genes
input_limmaFCViz <- 2 #FDR to filter genes
genomePlotly() # shows fold-changes on the genome
## Warning in eval(quote(list(...)), env): NAs introduced by coercion
## Warning in genomePlotly(): NAs introduced by coercion
input_nGenesBiclust <- 1000 #Top genes for biclustering
input_biclustMethod <- 'BCCC()' #Method: 'BCCC', 'QUBIC', 'runibic' ...
biclustering.out = biclustering() # run analysis
input_selectBicluster <- 1 #select a cluster
biclustHeatmap() # heatmap for selected cluster
input_selectGO4 <- 'GOBP' #Gene set category
# Read pathway data again
GeneSets.out <-readGeneSets( geneSetFile,
convertedData.out, input_selectGO4,input_selectOrg,
c(input_minSetSize, input_maxSetSize) )
results <- geneListBclustGO() #Enrichment analysis for k-Means clusters
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| adj.Pval | Genes | Pathways |
|---|---|---|
| 5.3e-93 | 260 | Response to abiotic stimulus |
| 1.1e-51 | 173 | Response to oxygen-containing compound |
| 4.1e-49 | 146 | Cellular response to stress |
| 5.1e-47 | 161 | Cellular response to chemical stimulus |
| 5.3e-46 | 185 | Response to organic substance |
| 3.6e-45 | 63 | Cellular response to hypoxia |
| 4.1e-45 | 63 | Cellular response to decreased oxygen levels |
| 4.1e-45 | 63 | Cellular response to oxygen levels |
| 1.1e-41 | 63 | Response to hypoxia |
| 2.8e-41 | 63 | Response to decreased oxygen levels |
input_mySoftPower <- 5 #SoftPower to cutoff
input_nGenesNetwork <- 1000 #Number of top genes
input_minModuleSize <- 20 #Module size minimum
wgcna.out = wgcna() # run WGCNA
## Warning: executing %dopar% sequentially: no parallel backend registered
## Power SFT.R.sq slope truncated.R.sq mean.k. median.k. max.k.
## 1 1 0.720000 1.5400 0.872 373.00 385.00 521.0
## 2 2 0.219000 0.4450 0.727 201.00 194.00 352.0
## 3 3 0.000261 -0.0124 0.692 125.00 113.00 262.0
## 4 4 0.192000 -0.3440 0.729 85.10 73.20 204.0
## 5 5 0.425000 -0.6040 0.808 60.70 51.00 164.0
## 6 6 0.596000 -0.7820 0.877 44.90 37.00 134.0
## 7 7 0.688000 -0.9280 0.925 34.10 27.50 111.0
## 8 8 0.736000 -1.0400 0.940 26.40 20.70 93.5
## 9 9 0.754000 -1.1500 0.951 20.80 15.90 79.3
## 10 10 0.769000 -1.2100 0.967 16.60 12.30 67.8
## 11 12 0.810000 -1.3200 0.975 11.00 7.69 50.5
## 12 14 0.805000 -1.4100 0.975 7.52 5.03 38.4
## 13 16 0.813000 -1.4700 0.963 5.29 3.22 29.6
## 14 18 0.826000 -1.4900 0.981 3.82 2.16 23.2
## 15 20 0.813000 -1.5100 0.961 2.81 1.51 18.3
## TOM calculation: adjacency..
## ..will not use multithreading.
## Fraction of slow calculations: 0.000000
## ..connectivity..
## ..matrix multiplication (system BLAS)..
## ..normalization..
## ..done.
softPower() # soft power curve
modulePlot() # plot modules
listWGCNA.Modules.out = listWGCNA.Modules() #modules
input_selectGO5 <- 'GOBP' #Gene set category
# Read pathway data again
GeneSets.out <-readGeneSets( geneSetFile,
convertedData.out, input_selectGO5,input_selectOrg,
c(input_minSetSize, input_maxSetSize) )
input_selectWGCNA.Module <- 'Entire network' #Select a module
input_topGenesNetwork <- 10 #SoftPower to cutoff
input_edgeThreshold <- 0.4 #Number of top genes
moduleNetwork() # show network of top genes in selected module
## softConnectivity: FYI: connecitivty of genes with less than 15 valid samples will be returned as NA.
## ..calculating connectivities..
input_removeRedudantSets <- TRUE #Remove redundant gene sets
results <- networkModuleGO() #Enrichment analysis of selected module
results$adj.Pval <- format( results$adj.Pval,digits=3 )
kable( results, row.names=FALSE) %>%
kable_styling(bootstrap_options = c("striped", "hover")) %>%
scroll_box(width = "100%")
| adj.Pval | Genes | Pathways |
|---|---|---|
| 5.3e-93 | 260 | Response to abiotic stimulus |
| 1.1e-51 | 173 | Response to oxygen-containing compound |
| 4.1e-49 | 146 | Cellular response to stress |
| 5.1e-47 | 161 | Cellular response to chemical stimulus |
| 5.3e-46 | 185 | Response to organic substance |
| 3.6e-45 | 63 | Cellular response to hypoxia |
| 4.1e-45 | 63 | Cellular response to decreased oxygen levels |
| 4.1e-45 | 63 | Cellular response to oxygen levels |
| 1.1e-41 | 63 | Response to hypoxia |
| 2.8e-41 | 63 | Response to decreased oxygen levels |